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Abstract 

We investigate properties of a solar wind-like plasma including a secondary alpha particle popula¬ 
tion exhibiting a parallel temperature anisotropy with respect to the background magnetic field, using 
linear and quasi-linear predictions and by means of one-dimensional hybrid simulations. We show 
that anisotropic alpha particles can drive a parallel fire hose instability analogous to that generated 
by protons, but that, remarkably, the instability can be triggered also when the parallel plasma beta 
of alpha particles is below unity. The wave activity generated by the alpha anisotropy affects the 
evolution of the more abundant protons, leading to their anisotropic heating. When both ion species 
have sufficient parallel anisotropies both of them can drive the instability, and we observe generation 
of two distinct peaks in the spectra of the fluctuations, with longer wavelengths associated to alphas 
and shorter ones to protons. If a non-zero relative drift is present, the unstable modes propagate 
preferentially in the direction of the drift associated with the unstable species. The generated waves 
scatter particles and reduce their temperature anisotropy to marginally stable state, and, moreover, 
they significantly reduce the relative drift between the two ion populations. The coexistence of modes 
excited by both species leads to saturation of the plasma in distinct regions of the beta/anisotropy 
parameter space for protons and alpha particles, in good agreement with in situ solar wind observa¬ 
tions. Our results confirm that fire hose instabilities are likely at work in the solar wind and limit the 
anisotropy of different ion species in the plasma. 

Subject headings: Solar wind; plasma; instabilities; waves; numerical simulation 


1. INTRODUCTION 

Temperature anisotropies with respect to the direction 
of the local magnetic field are ubiquitous in the weakly 
collisional heliospheric plasma; particle distribution func¬ 
tions measured in situ by spacecraft at various heliocen¬ 
tric distances in the solar wind are not at thermal equi¬ 
librium and de partures from Maxwe ll ian are frequently 
observed fe.g.. iMarsch et al.l Il982bt iNeueebauer et al.l 
Il996t iGarv et ai.ll2006l iMatteini et al.ll2013af) . The en- 
ergy transport in non-thermal plasmas is an impor¬ 
tant topic for space and astrophysical environments, 
where microphysical processes generated by tempera¬ 
ture anisotropies can importantly affe ct and change 
the m acroscopic evolution of the system (IMatteini et al.l 
l2012f) . Instabilities driven by a temperature anisotropy 
are thought to be an important mechanism in weakly 
collisional or collisionless plasma, where Coulomb col¬ 
lisions are not able to maintain the thermodynamical 
equilibrium of the gas and other processes take place to 
regulate the evolution of the particle distribution func¬ 
tions. Direct mea surements of the ion properties in 
the s olar wind (e.g.. lHellinger et al.ll2006HMatteini et al.l 
I2007f) and in the magnetosphere (fSamsonov et al.l 1200711 
suggest that kinetic instabilities driven by an ion tem¬ 
perature anisotropy play a role in constraining the ra¬ 
tio of parallel and perpendicular temperatures, defined 
with respect to the local mean magnetic field. More¬ 
over, observations of electromagnetic spectra in the so¬ 
lar wind show some possible signatures of the presence 
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of fluctuations associated to those processes ([Bale et al.l 
120091 IWicks et al.l 1201^ . although their interpr etation 
needs some caution (cfr. IHellinger fc: Travnicekl 1201411 . 
Kinetic instabilities like mirror, ion-cyclotron, and fire 
hose are also likely to b e active in astrophys i cal systems, 
e.g. accretion disks (iSharma et al.l 1 20061: iKunz et al.l 


120141: iRiauelme et al.ll2015l: ISironi fc Naravan 2015 ) and 

the hot plasma in galaxy clusters ([Schekochihin et al.l 
120051: IMogavero fc Schekochihii] 12014) . an d in plasmas 
wher e magnetic reconn e ction i s at work (iDrake et ahl 
20101: iS choeflle r et al.l 120121: IMatteini et al.1 l2013bl : 


Gingell et al.ll2015l)~ 


Among the different instabilities, fire hose is expected 
to be the dominant process in regulating the temperature 
anisotropy of a weakly collisional plasma when the paral¬ 
lel temperature is larger than the perpendicular temper¬ 
ature. Note that this condition is naturally obtained in 
systems where the magnetic field is spatially modulated 
and decreasing with distance (as for example expand¬ 
ing plasmas, like the solar and other magnetized stel- 
lar winds) due to part icle magnetic moment conservation 
(IMatteini et al.|[2012l) or where field aligned beams - con¬ 
tributing to an excess of the total parallel temperature - 
are gen erated, as typically observed in reconnec tion ex¬ 
hausts (iHoshino et al.|[T998l: iGosling et ani2005L Hietala 
et ah, submitted). 

Previous studies of fire hose instabilities going beyond 
the fl uid approximation have focussed on both electrons 
(e.g.jiLi fc: Habba]|l200CII: lMessmeril200^ iPaesold fc Benz 
2003; iGamporeale fc BurgessI 20081: Hellinger et al 


2014) and pro t ons (e.g.. iQuest fc Shapirol 19961: 
Gary et ^ |199^ [Hellinger e t al. 120031 : IMatteini et al l 


20061: IHellinger fc TravniCekI 1200^ . These studies are 


supported by in situ observations in the solar wind. 
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for protons (jKasper et al.l l2002t iHellinger et al.l I2006D 
and electrons (jStverak et al.l I2nn8jl . However, much 
less attention has been addressed, to our knowledge, 
to the case when the fire hose instability is driven by 
other species of the plasma. While the alpha particle 
abundance in the solar wind is only few percent, they are 
heavier and usually both faster and hotter than protons 
so that they carry important fractions of the solar wind 
momentum and energy. Their deceleration with respect 
to the protons may serve as an important source of 
energy for the proton heating (ISc hwartz fc Mysc hlllQSSl 
IHellinger fc Travmcekil2013HVerscharen et al.ll201^ 

Recent observations of t he solar wind alpha parti¬ 
cle distribution functions (IMaruca et al.l I2012I1 . have 
shown that temperature anisotropies of alpha particles 
also appear to be constrained by boundaries in plasma 
beta/anisotropy space, suggesting the role of micro- 
insta bilities in the regulation o f alpha properties. More¬ 
over, iBourouaine et all (1201311 have shown that alpha- 
protons drifts may also be influenced by these processes. 

How this happens and how different instabilities driven 
by different species simultaneously act in a plasma and 
stabilize particles to different marginally stable states has 
not been investigated in detail yet and it is the aim of 
this paper to do it by means of both linear and quasi- 
linear theoretical calculations and fully nonlinear numer¬ 
ical simulations. We show results from 1-D numerical 
hybrid simulations that describe the development of the 
parallel fire hose, which is the dominant instability when 
a composition compatible with the solar wind plasma is 
adopted (see below), in the presence of anisotropic al¬ 
phas with variable properties. Nevertheless, results ob¬ 
tained in this particular framework may be applied, with 
a suitable change in the parameters adopted, to any other 
space or astrophysical plasma. 

The paper is organized as follows: in Section [5] we 
present linear theory predictions of the parallel fire hose, 
focussing in particular on the differences between proton 
and alpha contributions in driving the instability and the 
influence of a velocity drift between the species; in Sec¬ 
tion [3] we discuss the feedback of the instability on the 
plasma, deriving the ion parallel and perpendicular heat¬ 
ing rates and momentum exchange expected on the basis 
of a quasilinear approach; in Section U] we show results 
from hybrid simulations for different parameter combi¬ 
nations, including both the case of purely alpha driven 
fire hose as well as cases with mixed ion contributions. 
Finally in Section [S] we summarize and discuss the results 
obtained, comparing linear and quasi-linear predictions 
with the fully nonlinear evolution found in simulations 
and in Section [5] we present our conclusion, highlighting 
implications of our study for solar wind and astrophysical 
plasmas. 

2. LINEAR THEORY 

We start our analysis with linear predictions, obtained 
solving the dispersion of a fully kinetic Vlasov system. 
The dispersion and stability of the plasma depends on 
many parameters and the presence of three popula¬ 
tions complicates further the picture with respect to the 
case of a pure electron-proton plasma. For a detailed 
linear theory investigation of the role of alpha parti¬ 
cles a nd their i nfluen ce on the dispersion see for ex¬ 
ample iMarucal (|2012(1 . Two Are hose instabilities are 



Figure 1. Fire hose unstable region in the presence of alpha parti¬ 
cle temperature anisotropy; contours of the 7m = 10“^ maximum 
growth rate level are shown for a case with an alpha-proton drift 
Vap = 0.6u^ (solid line), and with v^p = 0 (dashed line). The pro¬ 
ton driven (ria = 0) and fluid CGL fire hose thresholds are shown 
as references in dash-dotted and dotted lines, respectively. 

in general active in plasmas with T± < T|| (see Ap- 
pendix for symbol definitions): the paral lel Are hose (e.g., 
IQuest fc ShapirolllQQGtlGarv et alli998ll . with maximum 
growt h rate at parallel propagation, and t he oblique Are 
hose ([Hellinger fc Matsumotal2000L l2001h . which has a 
maximum for oblique wave-vectors. In a plasma with 
proton anisotropy the parallel instability is typically 
dominant for fHiip < 10. The same behavior is predicted 
by the linear theory when the instability is driven by 
temperature anisotropy in alpha particles, with a prefer¬ 
ential excitation of fluctuations parallel to the magnetic 
field. For this reason, in this work we specialize to the 
case of the dominant parallel Are hose instability. 

Figure [T] shows the Are hose unstable region in the pa¬ 
rameter space (/3||s) Tj-s/r'iis) derived from Vlasov linear 
theory for diflerent plasma conflgurations. Contours dis¬ 
play the 7m = 10“^ maximum growth rate level, which 
can be reasonably taken as an estimation of the threshold 
(or rather a marginal stability condition) of the instabil¬ 
ity (e.g.. iMatteini et al.ll2n0ffl . Solid and dashed lines 
are associated with anisotropic bi-Maxwellian alpha par¬ 
ticles, with density = O.Ohrie; for these computations 
(where the axis refers to alpha particles) protons and 
electrons are isotropic, and /3p = /3e = 1. The solid line 
shows the maximum growth rate of the alpha Are hose 
instability when an alpha-proton drift Vap = O.Gva is as¬ 
sumed in the plasma; not e that this is a typical v alue in 
the fast solar wind (e.s.. iNeugebauer et al.lflM^ . The 
dashed line refers to the same level for the case with 
Vap = 0. As a reference, the dotted line shows the CGL- 
fluid Are hose threshold: fiu —13± = 2. Note flrst that the 
unstable region predicted for anisotropic alphas extends 
to much smaller betas with respect to the fluid case. This 
remains true also when smaller values are assumed for the 
proton and electron betas, which contribute to the total 
plasma pressures, underlining that anisotropic alpha par¬ 
ticles can drive a kinetic Are hose instability also when 
the total plasma beta is <C 2, thus entirely stable with re¬ 
spect to the corresponding fluid instability. Such a linear 
prediction is conflrmed by numerical simulations, as we 
will show in what follows. For the choice of parameters 
used here, the alpha-driven unstable region is also larger 
than the analogous proton parallel Are hose instability 
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Figure 2. Instability growth rates for /3p = 3 and /3a = 2 and anisotropic protons Ap = 0.5 and alphas Ac = 0.4. Left panel refers 
to forward propagating waves and right panel to backward propagating modes. The red line shows case with Vap = 0, while dashed line 
refers to different values of the alpha-proton drift Vap = 0.3, 0.6, O.Qnyi. The dash-dotted line shows the growth rate for the same alpha 
parameters but isotropic protons. 


region predicted for typical solar wind conditions (see 
(IMatteini et al.l[2006f) l when alphas are not taken into 
account, which is approximatively bounded by ,5||p > 1 
and is shown as the dash-dotted contour. 

Second, the unstable region is wider in the case of 
a non-zero alpha-proton relative drift, and in general 
the growth rate significantly depends on the value of 
Vap. Note that the presence of a relative drift between 
species also influences the spectrum of waves excited by 
the instability in terms of their direction of propagation, 
and then modulates the associate d scattering of parti¬ 
cles (e.g.. iVerscharen et al.|[^13a|l . In the absence of a 
velocity separation between protons and alphas a balance 
in the generation of forward and backward propagating 
modes is expected. By contrast, in case of an alpha- 
proton relative drift, this symmetry is broken; if Vap ^ 0 
then dominant unstable waves of each species are ex¬ 
pected to propagate in opposite directions. Resonant fire 
hose is driven by the anomalous-doppler resonance (e.g. 
IKennel fc Scarilll968h : 

uj - k\ivi\ =-fls , ( 1 ) 

where w and are respectively the frequency and the 
parallel wave-vector of the unstable wave, W[| is the 
parallel velocity of the resonant particles, and fig is 
the cyclotron frequency of species s (Note that in this 
resonance the Doppler-shifted frequency has opposite 
sign with respect to the standard cyclotron resonance: 
oj — k\\v\\ = Dg). Resonant parallel fire hose instability 
destabilizes right-handed magnetosonic/whistler waves, 
and since particles satisfy the resonant condition m with 
kv\\ > 0, we expect then the unstable modes to be pref¬ 
erentially generated in the direction of the drift of the 
corresponding species. For alphas drifting ahead protons 
along the magnetic field, this means they support mostly 
forward propagating waves, while protons tend to excite 
more backward propagating modes. 

This scenario is confirmed by linear theory and it can 
be seen in Fig. where the two panels show predictions 
for unstable forward (left panel) and backward (right 
panel) propagating modes. The red line (same in both 
panels) shows the predicted growth rates in case when 
no drift is present and both species are anisotropic and 
unstable {Ap = 0.5 and Aq,= 0.4; /3p = 3 and /3a = 2); in 
this case symmetry of the system produces equal amount 
of waves in both direction. Also note that the presence 
of anisotropic protons significantly enhances the growth 


rate of the instability with respect to the case when only 
alphas are anisotropic, shown in dash-dotted line, and 
that this is true also in cases when protons are anisotropic 
but stable with respect to the fire hose. It is interesting to 
observe that, for the parameters used in Figure HJ there 
are two distinct peaks in the growth rate as a function 
of wave-vectors. These correspond to distinct fire hose 
modes driven by proton and alpha anisotropies, with the 
lower k peak associated with the alphas and the other 
with the protons. This is a general feature characteris¬ 
tic of the regime investigated in this work, and is also 
confirmed by numerical simulations. 

As expected when a relative drift is introduced 
{vap/vA = 0.3, 0.6, and 0.9, shown by dashed lines), the 
symmetry between left and right panels is broken. As Vap 
is increased the peak corresponding to alpha fire hose is 
enhanced for forward waves (same direction as the drift of 
alphas in the electron/plasma frame), while it is slightly 
weakened for backward waves. The opposite is observed 
for protons, for which the corresponding peak increases 
for backward waves and decreases for forward, owing to 
their negative drift in plasma frame. 

3. QUASILINEAR PREDICTIONS 

A further step in the analysis can be done estimat¬ 
ing quasiline ar heating rate s assoc iated to the linear 
predictions. iHellinger et al.l (l2013f) estimated the ini¬ 
tial quasilinear parallel and perpendicular heating rates 
(per mode) by applying the quasilinear diffusion opera¬ 
tor (containing only the contribution of the given mode) 
on the initial bi-Maxwellian particle velocity distribu¬ 
tion functions and calculating the corresponding mo¬ 
ments. Following this approach, we have then calculated 
the quasilinear plasma response to fluctuations gener¬ 
ated by unstable ions, for parameters studied in Fig. ^ 
Aa = 0.4, /3p = 3, Pa = 2, and as a function of the range 
of unstable modes k. Fig.[3]shows the real and imaginary 
part of the complex frequency from linear theory (top), 
and the proton (middle) and alpha (bottom) quasilinear 
predictions. In the second and third panels, the dashed 
line shows the perpendicular heating rate, the solid line 
the parallel heating rate (negative values mean cooling), 
the blue line shows the total kinetic (thermal) energy 
and the red dotted line shows the total energy of the sys¬ 
tem (kinetic-kmagnetic) that is conserved to a very good 
approximation. Heating rates are normalized to Eem'-Op, 
where E^m is the electromagnetic energy of a single in- 














4 


Matteini et al. 








Figure 3. Linear and quasi-linear predictions for the same plasma parameters as in Fig. [2] Aq, = 0.4,/3p = 3, /3 q; = 2, and with Ap = 1 
(left panels, dash-dotted line case in Fig. and with Ap = 0.5 (right panels, red line case in Fig. IS. Top: real frequency (dot-dashed) 
and growth rate (red) as a function of the unstable wave-vectors k. Quasilinear predictions of heating rates dEjdt as a function of k for 
protons (middle) and alphas (bottom) are also shown; dashed and solid lines refers to parallel and perpendicular directions, respectively. 
The blue line shows t he total heating rate and the dotted red line the total energy of the system. Heating rates are normalized following 
IHellinger et al.l <120131) and negative values correspond to particle cooling. 


teracting mode fsee IHellinger et HI (j201,SD for details). 

A case when only alpha particles are anisotropic and 
drive the instability - corresponding to the dash-dotted 
line in Fig. {Ap = 1) - is shown in the left panels of 
Fig. [3l This calculation results in a perpendicular heat¬ 
ing and a parallel cooling for alphas (middle panel), as 
expected for fire hose instability; such an energy transfer 
leads to the isotropization of the initial T± < T|| condi¬ 
tion, owing the self-regulating behavior of the instabil¬ 
ity. It is interesting to note that the same response is 
predicted also for protons (bottom-left panel); quasilin¬ 
ear approach predicts perpendicular heating and parallel 
cooling for the isotropic and stable protons, as a response 
to the generation of fluctuations by the alphas. As it will 
be discussed in the next session, this prediction is con¬ 
firmed by fully nonlinear numerical simulations. This 
also means that part of the free energy extracted from 
alpha particles is transferred to protons, via diffusion by 
the unstable fire hose modes. Note that this also leads 
to the generation of a proton temperature anisotropy 
T± > T||. 

The pictures is further complicated when both species 
are initially anisotropic and contribute to the instability. 
As discussed previously, linear theory predicts distinct 
unstable waves associated with the two populations. In 
this case particles can also resonate with modes driven 
from the other species. This is shown in the right panels 
of Fig. [3] for the same parameter combination as the red 
line in Fig. [3] {Ap = 0.5 and Vap = 0). By comparing 
the three panels, it can be noticed that each species is 
isotropised mostly in correspondence to the peak of the 
growth rate associated with their anisotropy. However, 
protons are significantly anisotropically heated/cooled 


also by fluctuations that are driven by alphas, at lower 
k. By contrast, modes generated by protons {k > 0.3) 
have little influence on alphas. 

If a drift velocity is present, then also a slowing down 
of the alpha-proton drift is expected, together with the 
change in the thermal energetics described above. The 
rate of the change in the kinetic energy associated with 
the proton-alpha drift induced by fire hose can be anal¬ 
ogously derived from the quasilinear theory (Hellinger 
et al. 2013). In this case, the quasilinear calculation 
predicts a decrease of the drift between the two species, 
with a fraction of the kinetic energy associated with the 
differential motion converted into thermal energy (and 
also into wave generation), leading to further anisotropic 
heating of both species. These theoretical quasilinear 
predictions, however, get complicated when too many pa¬ 
rameters are taken into account, and, moreover, they are 
only valid initially when the particle distribution func¬ 
tions are close to bi-Maxwellian. For a more complete 
description of the plasma evolution, numerical simula¬ 
tions are a more suitable tool. 

4. SIMULATIONS 

We have performed one-dime nsional (l-D) nu merical 
simulations using a hybrid code (lMatthewslll994H which 
treats protons as particles and electrons as a massless 
charge neutralizing fluid with a constant temperature. 
In the code units of space and time are the ion inertial 
length and the inverse proton cyclotron frequency, re¬ 
spectively. See the Appendix for definition of code units. 
We use a simulation grid of 1200 collocation points, with 
spatial resolution Ax = 1 and 10^ particles per cell for 
each species. 

We have performed several runs, with different initial 
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Figure 4. Example of alpha fire hose instability (Run C). Top 
left: alpha temperature anisotropy evolution; Top right; magnetic 
energy evolution; Bottom left; spectrum of fluctuations at t = 
150 and linear theory prediction of instability growth rates (red); 
Bottom right: alpha particle distribution at t = 0 (red line) and 
t = 150. 

conditions corresponding to the different cases discussed 
through the previous sessions and summarized in Ta¬ 
ble m Run A-C focus on the case when only alpha parti¬ 
cles are anisotropic over different beta regimes, Run D-F 
describe configurations where both protons and alphas 
are anisotropic, in order to study the coexistence of the 
modes driven by both ions, and finally Run G-I include 
a relative drift between the species. The main results of 
these simulations are presented and discussed over the 
next sections. 


4.1. Isotropic protons 

We start with a configuration where only alpha par¬ 
ticles are initially anisotropic. Figured provides an 
overview of the dynamics observed for Run C, which has 
plasma parameters c onsistent with solar wind conditions 
(|Mariica et al.lf2(n^ . Top panels show the temporal evo¬ 
lution of the anisotropy (left) and of the magnetic fluctu¬ 
ations (right). A linear phase of the instability is visible 
in the interval 50 < f < 100, as suggested by the expo¬ 
nential growth of the magnetic fluctuating energy, which 
then saturates after t ^ 120. As the unstable fluctu¬ 
ations develop, they start scattering particles reducing 
their initial anisotropy. At saturation the plasma has 
reached a marginally stable state with T±a/T\ia ^ 0-4. 
The bottom left panel shows the spectrum of the fire 


Table 1 

Initial parameters for Hybrid 1-D simulations 


Run 

P\\p 

Ap 

P\\a 

Acx 

'^ap / '^A 


0\\tot 

A 

0.1 

1 

0.4 

0.1 

0 

0.01 

0.41 

B 

2 

1 

2 

0.2 

0 

1 

5 

C 

6 

1 

2 

0.2 

0 

3 

11 

D 

1.5 

0.2 

0.7 

0.15 

0 

0.5 

2.7 

E 

2 

0.3 

1 

0.25 

0 

0.5 

3.5 

F 

3 

0.5 

2 

0.4 

0 

1 

6 

G 

1 

1 

0.5 

0.5 

0.9 

0.5 

2 

H 

2.5 

0.5 

0.6 

0.5 

0.9 

0.5 

3.6 

I 

3 

0.5 

2 

0.4 

0.9 

1 
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Figure 5. (Top) Evolution of parallel and perpendicular proton 
temperatures in Run C; (bottom) evolution of the power of mag¬ 
netic fluctuations as in Fig. |4] 


hose fluctuations at t = 120, at saturation. The red 
line corresponds to the growth rate predicted by the lin¬ 
ear theory. There is a very good agreement between the 
range of unstable modes predicted by the theory and the 
spectrum observed in the simulation after saturation. Fi¬ 
nally, in the right bottom panel the alpha particle distri¬ 
bution function after saturation is shown; the shape of 
the distribution at f = 0 is also shown as a red contour. 
As a consequence of the resonant wave-particle interac¬ 
tion between alpha particles and fire hose fluctuations, 
pitch-angle scattering signatures in velocity space are ob¬ 
served producing a significant deviation from the initial 
bi-Maxwellian. The deformation of the distribution is 
similar to t he one observed for p rotons in the case of 
the parallel (IMat teini et al.ll2006D and oblique fire hose 
(|Hellinger fc Travnicekl I2008lh Note that, as discussed, 
the position of the unstable spectrum in k-space and the 
signatures in the particle distributions depend on reso¬ 
nances and are then a function of the specific plasma 
parameters taken as initial conditions; however the same 
qualitative global behavior of the instability shown in 
Fig. |3]is found also in simulations with different (lower) 
initial alpha and proton betas (Run A-B). Note that in 
Run A the total plasma beta is smaller than 1, corre¬ 
sponding then to a condition theoretically stable with 
respect the CGL-ffuid fire hose (see Figure[T]). Although 
plasma parameters in this case are not realistic in terms 
of solar wind conditions, resulting in an excessively large 
alpha-proton temperature ratio, more generally Run A 
highlights the kinetic nature of the instability. 

It is interesting to investigate the feedback on the back¬ 
ground dominant proton population from the dynamics 
driven by the secondary alpha population. As shown in 
the previous section, the quasilinear approximation pre¬ 
dicts an anisotropic heating/cooling of the proton distri¬ 
bution by the modes excited by the alphas; in particular 
the generation of a proton temperature anisotropy is ex¬ 
pected. This is confirmed by Figure (5] where the solid 
and dashed lines in the top panel refer to the evolution 
of the parallel and perpendicular proton temperatures, 
respectively. The development of a proton anisotropy is 
very well correlated to the growth of unstable fluctua¬ 
tions driven by the alpha fire hose instability (bottom 
panel). Moreover, the temperature evolution predicted 
by the quasilinear calculation (perpendicular heating and 
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parallel cooling) describes well the evolution of the tem¬ 
perature profiles during the initial phase of the instabil¬ 
ity. After saturation {t ^ 150), well beyond the validity 
of the linear and quasilinear approximations, a parallel 
proton heating is observed, which progressively reduces 
the proton anisotropy developed during the linear phase 
of the instability; this is likely related to the late damping 
of the fire hose fluctuations, as a slight decrease in the 
magnetic field power after saturation is observed (bottom 
panel). 

4.2. Alphas and protons unstable 

We consider now the situation when both species are 
initially anisotropic, which constitutes a more realistic 
description of the solar wind regime and possibly of 
weakly collisional astrophysical plasmas in general. In 
this configuration, when starting with both species pre¬ 
dicted unstable by the linear theory, we expect to observe 
a co-existence of the proton and alpha fire hose instabil¬ 
ities. Figure [51 top panel, shows the spectrum of fluc¬ 
tuations for Run F, for which a double-peaked spectrum 
of fluctuations, as in Fig. [21 is expected. The spectrum 
observed in the simulation is in good agreement with the 
linear prediction of the unstable modes (red line), show¬ 
ing distinct contribution from proton and alpha unstable 
modes. The blue line shows a smoothed profile of the 
spectrum in the simulation, and highlights the presence 
of a double peaked profile. As discussed in section |2l the 
generation of unstable fluctuations at larger wavelength 
is associated with the modes driven by the alphas (con¬ 
sistent with the lower k single peaked spectrum observed 
in Fig. m when protons are stable), while smaller wave¬ 
lengths correspond to proton modes. Moreover, distri¬ 
butions are symmetric with respect to the magnetic field 
direction, so that modes propagating in both directions 
are developed in the simulation, as expected; the bot¬ 
tom panel of the figure shows a profile of the magnetic 
field component Bz as a function of time, highlighting 
the presence of counter-propagating wavefronts approxi- 
matively propagating at -^va (white trajectories in solid 
line), driven by the instabilities after t ^ 40. 

In order to investigate further the evolution of the in¬ 
stability when the two species are both unstable, it is use¬ 
ful to use the (/Sp, Tj_s/T||s) plane. FigurejTjshows trajec¬ 
tories in this parameter space for three simulations with 
different initial conditions (Run D-E-F: solid, dashed and 
dotted line, respectively). Note that parameters of these 
runs, chosen in order to emphasize the simultaneous evo¬ 
lution of both species close to fire hose instability thresh¬ 
olds, lead to an initial alpha-proton temperature ratio 
which is larger than typical solar wind conditions. How¬ 
ever, we have verified that this choice does not change 
significantly the results and does not affect the qualita¬ 
tive evolution of the plasma with respect to simulations 
with a more realistic temperature ratio. The solid con¬ 
tours correspond to the alpha fire hose unstable region, 
as in Figured! while the dashed contours show the insta¬ 
bility threshold for the instability driven by protons. The 
red and blue lines indicate the evolution path of alphas 
and protons, respectively, and symbols show the initial 
conditions. Both species are initially unstable and then 
become more isotropic during the simulations; this cor¬ 
responds to a migration towards the top-left corner in 
the parameter space. Remarkably, each species is sta- 




t 


Figure 6. Run F; (top) Spectrum of the unstable modes at t = 80. 
Red line shows the linear prediction for the same parameters as in 
the run (See Fig. o, while the blue line is a smoothed profile of 
the simulation data, (bottom) Profile of the Bz component of the 
magnetic field as a function of time. The white straight lines show 
the propagation of forward and backward waves with phase speeds 
±VA- 

bilized according to its own marginal stability, thus in 
different regions of the parameter space: protons satu¬ 
rate at /3||p > 1, while alphas are pushed towards their 
threshold (solid contours) that extends also to /3||„ < 1. 
The final stable state found after the nonlinear stage in 
the simulation is in reasonable agreement with marginal 
stability given by linear theory. Note however that due 
to the resonant interaction between particles and unsta¬ 
ble modes, distributions may deviate sig nificantly from 
Maxw ellian at saturation; as discussed in iMatteini et al.l 
(120061) , this can lead to the fact that the plasma saturates 
in the parameter space before reaching the theoretical 
linear thresholds that are computed for bi-Maxwellian 
distributions. It is then not surprising that the agree¬ 
ment between simulations and linear theory in Fig. |7|is 
not exact, although it captures the global trend of the 
stabilization. 

4.3. Role of alpha-proton drift 

We finally consider the role of an alpha-proton drift 
in the dynamics of the instability. As shown in section 
m the presence of a relative streaming between species 
introduces an asymmetry in the system, owing to each 
species’ drift with respect to the plasma (electron fluid) 
frame. As a consequence, we expect to observe in the 
simulation a preferred direction of propagation in the 
modes driven by each species, following the direction of 
their drift in the plasma frame. To confirm this scenario, 
we present in this section results from runs where an 
initial drift of the order of the Alfven speed is taken into 
account. This i s a characteristic oft e n observed in the 
fast solar wind ( iMarsch et al.|[l982at iNeueebauer et aH 
Il996t [Goldstein et al.ll200(^ 

Figure [5] shows the results for Run I, which has the 
same initial conditions as Run F, but with Vap = 0.9?;^. 
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Figure 7. Evolution of anisotropies for alphas (red) and proton 
(blue); symbols show the unstable initial conditions for 3 different 
simulations: Run D (diamonds), Run E (stars), and Run F (trian¬ 
gles) . The dashed and solid contour lines show the constant growth 
rates levels of proton and alpha parallel fire hoses predicted by the 
linear theory, respectively. 

The top panel shows the evolution of temperature 
anisotropies; protons are affected first, at t ^ 30, lead¬ 
ing to a small change in their anisotropy, but sufficient 
to stabilize them. The activity driven by the protons 
has little influence on the alpha anisotropy, consistent 
with the quasi-linear prediction of Fig. [3l Shortly af¬ 
ter, from t ^ 50-60, alpha particles are isotropized more 
significantly, corresponding to the phase when fluctua¬ 
tions associated with the alphas start to play a role. A 
second change in the proton anisotropy (solid line), is 
also observed at t ^ 80, during the isotropization of al¬ 
phas. This phase is analogous to the anisotropic proton 
heating seen in Figure [5] (only alphas unstable) and it is 
also consistent with the quasilinear prediction discussed 
in section [3] and shown in Figure |3] 

During the activity of the instability the alpha-proton 
relative drift is also influenced. The red dash-dotted line 
in the top panel shows the evolution of Vap as a function 
of time. A first weak decrease can be observed when pro¬ 
tons become unstable, while the largest variation occurs 
later, when the alpha fire hose instability takes place; at 
the end of the run, the drift between the two ion popu¬ 
lation is reduced from Vap = O.Oua to Vap = O.Gva- We 
have checked that this behavior, which is in agreement 
with the quasilinear predictions, is confirmed by simu¬ 
lations where protons are isotropic (Run G) and with 
< 1 (Run H). 

The lower panel of Fig. [5]shows the profile of the trans¬ 
verse magnetic field component Bz{x), inside the region 
X = [0,400] and as a function of time. The white solid 
lines correspond to the trajectory of waves propagating 
with phase velocity Tua- Remember that, according to 
Figure the linear theory predicts for protons a domi¬ 
nance of backward propagating fluctuations, while alphas 
are expected to excite mainly forward propagating waves. 
Consistent with the expectations and the evolution of the 
anisotropies shown in the top panel, we observe the pres¬ 
ence of backward waves during the first phase of the run, 
when the instability is firstly driven by protons, due to 
their initial negative drift. After f ~ 70, when the macro¬ 
scopic dynamics begins to be controlled by the activity 
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Figure 8. Fire hose instability in the case of an alpha-proton 
drift Vap = 0.6da (Run I); in the top panel the solid line shows 
the alpha anisotropy, while dashed line refers to proton anisotropy, 
and red dot-dashed to the alpha-proton drift. In bottom panel, the 
profile of the transverse magnetic filed Bz {x) as a function of time 
is shown. The solid lines corresponds to the trajectory of wave 
propagating with ±ua. 

of alpha particles, forward modes emerge and dominate 
the magnetic field profile at later times. 

Finally, we consider the evolution of the ion total tem¬ 
peratures Tg during the instability. Fig.[3]shows the tem¬ 
poral profiles of Tp (top) and (bottom) normalized to 
the initial value, for Run F (dashed) and Run I (solid). 
When a relative drift between ion species is present, we 
observe a weaker cooling of the plasma during the fire 
hose activity. This is due to the fact that part of the ki¬ 
netic energy associated to the drift motion is converted 
into thermal energy (heating) of each species by the in¬ 
stability. This partially compensates the loss of particle 
energy produced by the instability and results in a heat¬ 
ing of each ion species. The effect is stronger for protons; 
moreover, during the nonlinear phase of the instability in 
Run I, after t ^ 120, we observe a further energization, 
leading to a net proton heating at expenses of the initial 
relative drift. 


5. DISCUSSION 

Properties of the alpha fire hose instability are found 
to be qualitatively similar to those of the instability gen¬ 
erated by the protons. Fig. |4] summarizes the global dy¬ 
namics of the instability: the exponential growth of un¬ 
stable modes extending over a finite range of k-vectors, 
the consequent reduction in anisotropy of particles by 
pitch-angle scattering and the associated deformation of 
the distribution function, developing characteristic non- 
thermal features. 

Differences with respect to the lighter protons are how¬ 
ever observed in the threshold of the instability and in 
the spectrum of unstable modes: fluctuations at larger 
wavelengths than in the case of protons are excited and 
the instability can be triggered at reasonably low betas 
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Figure 9. Temporal variation of total ion temperatures, normal¬ 
ized to initial values, during fire hose instability for a simulation 
with (Run I, solid) and without (Run F, dashed) proton-alpha drift. 
Top panel shows the evolution of the total temperature of protons, 
bottom panel of alphas. 

for typical solar wind and space plasma conditions (see 
Figure [T]). These properties can be qualitatively inter¬ 
preted in terms of the resonant interactions that drive the 
fire hose instabili t y in a kinetic regime. As discussed by 
iVerscharen et al.l (j2018aD . both protons and alpha par¬ 
ticles become fire hose unstable for Vth,s ^ va, where 
Vth,s is the thermal speed of each species. Due to the dif¬ 
ference in mass this translates into different conditions 
in terms of each species’ beta - for a given plasma com¬ 
position. Despite the fact that dispersion relations in a 
plasma with multi-species anisotropies and drifts is re¬ 
ally more complicated, this offers a simple argument to 
appreciate why in the case that fire hose is driven by al¬ 
phas, and for solar wind-like abundances, the parallel fire 
hose unstable region can extend to /3 ||q, < 1 with a sig¬ 
nificant departure from the fluid non-resonant threshold, 
as shown in Figure [1] Moreover, unlike the fluid non¬ 
resonant fire hose, the instability can be triggered also 
when the total beta of the plasma is well below unity 
(e.g.. Run A). Note that thi s can be relevant al so for as- 
trophysical and laboratory (jCarter et al.l 1201^ systems 
characterized by small beta regimes, where we expect 
heavier ions to be possibly more efficient in generating 
fire hose fluctuations, since protons are substantially fire 
hose stable for /3||p < 1. 

We have found that growth rates of the alpha fire hose 
are sensitive to several parameters; the presence of a 
proton anisotropy - even not necessarily large enough 
to drive the protons unstable - leads to a significant in¬ 
crease in the alpha instability growth rate (Figure (2), a 
trend w hich seems to be als o confirmed by in situ obser¬ 
vations (|Maruca et al.l[201^ . The behavior of the plasma 
is in particular modified by the presence of a relative 
drift between species. This breaks the symmetry in the 
propagation of unstable modes, leading to a net flux of 
forward/backward waves, depending on the drift of the 
species that drives the instability. This aspect is relevant 
for space plasmas, where differential streaming between 
the different components of the plasma is ubiquitously 
observed. As a consequence the parallel fire hose is ex¬ 
pected to produce a preferential generation of inward or 
outward waves depending on the species that drive the 
instability, thus offering a possible observational signa¬ 
ture for the identification of fluctuations generated by 
specific components of the plasma. 


On the other hand, we have shown that the fire hose 
instability plays also a role in regulating/reducing the 
differential speed itself. This is known in the case of 
fire hose driven by protons and has been proven pre¬ 
viously in numerical simul ations (iHellinger fc Travnicetd 
l2006t fMatteini et al.ll201^ . while we have found that in 
the case of an instability driven by the alphas, the effect 
can be even more significant. 

The parallel fire hose and the (parallel) magnetosonic 
beam instability can be considered as limiting cases of 
one instability driven by the same physical mechanism, 
the exce ss of kinetic/thermal ene rgy along the magnetic 
field (cf.. IVerscharen et al.ll2013b|l . Both instabilities lead 
to a significant reshaping of particle distributions. While 
beam-type instabilities with isotropic populations act to 
reduce the relative streaming between the ions down to 
the a value which is typically ~ va, the presence of tem¬ 
perature anisotropies can lead to a reduction o f the drift 
below that threshold (jVerscharen et al1l2013all . Consis¬ 
tently with that, in our simulations, as a result of the fire 
hose activity, the alpha-proton drift is reduced well below 
the Alfven speed {vap ^ 0.6va i n Figure Ell, in agreement 
with solar wind obs ervations (jNengebauer et ahl 119961: 
iGoldstein et al.lF2fllflf l. This further suggests that an in¬ 
stability driven by temperature anisotropy of alphas can 
play a role in regulating the alpha-proton drift in the so¬ 
lar wind. These results confirm that a complex interplay 
is possible between temperature anisotropy and relative 
drifts in space plasmas, and this needs more specific fu¬ 
ture investigations. 

6. SUMMARY AND CONCLUSION 

In summary, we have investigated the properties of the 
parallel fire hose instability generated by an alpha par¬ 
ticle pressure anisotropy. We have shown, for the first 
time, simulations of the linear and nonlinear phases of the 
instability, and we have compared the results with linear 
and quasilinear theory predictions. This study has led 
to several significant conclusions. First there is a good 
agreement in trends predicted by the linear and quasi¬ 
linear theories and the fully nonlinear results obtained 
during the early phase of the simulations. Quasilinear 
calculations correctly predict the evolution of the tem¬ 
peratures as a response to the fluctuations generated by 
the instability, leading to transverse heating and paral¬ 
lel cooling for all ions. Remarkably, the same behav¬ 
ior is observed also in the case when the main species 
(protons) is stable and isotropic, as a consequence of the 
activity driven by the less abundant alphas. This demon¬ 
strates that instabilities driven by secondary populations 
are also important for the global plasma energetics and 
that they can play a role in shaping the properties of the 
main ion population. Our simulations suggest that typi¬ 
cal solar wind abundances (^ 5%) are sufficient to trigger 
alpha fire hose growth rates that lead to dynamics with 
significant impact on the protons. 

It is important to highlight that the role of instabilities 
driven by alpha temperature aniso tropies is confirmed 
by direct solar wind observations (iMariica et al.l I2012H 
demo nstrating that, similarly to protons (|Hellinger et al.l 
l2006jl , the histogram of observational counts in the plane 
{P\\cn Tj_q,/T||q,) at lAU shows apparent empirical bound¬ 
aries that correspond well to the expected thresholds of 
kinetic instabilities in the same parameter space. More- 
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over, when comparing the proton and alpha particle ob¬ 
servations, the latter show a clear shift towards lower 
betas in their distribution (Figure 2 in iMaruca et al.l 
(j2012ll l. This scenario is consistent with the linear the¬ 
ory prediction (see Fig. [1]), where the alpha unstable re¬ 
gion extends to significantly smaller betas compared to 
both the case when the instability is driven by protons 
and the classical fluid threshold. As discussed in the 
previous section this can be explained in terms of the 
different effects of resonances that drive the instability 
for protons and for heavier alpha particles. Such a be¬ 
havior is very well confirmed by our simulations, where 
it is clearly observed that alphas and protons saturate 
and are constrained by fire hose instabilities in distinct 
regions of the parameter space (/3i|, Tl/T||), thus in very 
good agreement with the in situ observations. Moreover, 
we highlight the fact that the solar wind different ions are 
bounded simultaneously in distinct regions of the param¬ 
eter space, in good agreement with the distinct fire hose 
theoretical thresholds of both protons and alphas. This 
constitutes a further strong indication, together with pre¬ 
vious studies (e.g.. iKaspe r et al.l 1 20021: l^llinger et al. 
gOOliy Marsch et al.ll200(it Bale et al.ll2009tTMatteini et al. 

2018ajl . that this instability is at work in limiting the 
temperature anisotropy of ion species in the solar wind. 

Our simulations also show that as a result of the scat¬ 
tering by the unstable fluctuations, the fire hose insta¬ 
bility driven by alphas influences the evolution of the 
alpha-proton drift. We suggest that this mechanism can 
play a role in regulating ion-beams in the solar wind ex¬ 


pansion, since sub-Alfvenic drifts are often observed in 
the fast wind, a signature consistent with the results of 
our simulations. Clearly more work is needed to test the 
role of the alpha fire hose instability in the expanding 
solar wind and its coupling with othe r processes and mi¬ 
croinstabilities (iHellinger et al.ir2ni5ll . 

Finally, note that our findings are primarily addressed 
to heliospheric investigations, where a direct comparison 
between in situ data and simulations may be achieved; for 
this reason we have adopted abundances and properties 
of the plasma that are typical of the heliospheric plasma. 
However, several astrophysical environments, including 
galaxy clusters and accretion disks, are expected to be 
environments where pressure aniso tropies play a role 
in regulating the plasma energetics (|Schekochihin et ^ 
l2nn5HSharma et'ani20n6[ l . Despite some work about pro- 
ton driven instability having been rec ently carried out 
(iKnnz et al.ll20l4lRiaiielme et al.ll2015[ l. we suggest that 
the contribution from other minor species of the plasma 
may be important as well. 

This work was supported by UK Science and Tech¬ 
nology Facilities Council grant ST/K001051/1. The re¬ 
search leading to these results has received funding from 
the European Commission’s Seventh Framework Pro¬ 
gramme (FP7/2007-2013) under the grant agreement 
SHOCK (project number 284515). PH acknowledges 
GACR grant 15-10057S. 


APPENDIX 

APPENDIX: SYMBOLS AND UNITS 

Here T and || denote the directions with respect to the ambient magnetic field Bq, Bq = |Bo| its magnitude, and 5B 
denotes the magnitude of the fluctuating magnetic field; i;|| and v± are velocity components parallel and perpendicular 
to Bo, respectively; t denotes the time. Subscripts s refer to different plasma species; p stands for protons, a for alpha 
particles. T±s and Tp are the perpendicular and parallel temperatures, respectively, and Ag is the species temperature 
anisotropy Tls/Ths- Vth,s = (ksTs/msY^'^ is the thermal velocity and /3p = 87rnsfcBP[|s/^o parallel plasma 

beta for a given plasma specie s of density Ug (fcs is the Boltzmann constant). 

Through the paper the adopted units of space and time are the ion inertial length c/uip and the inverse proton 
cyclotron frequency respectively, where flp = eBo/mpC and Wp = (47rne^/mp)^/^ is the proton plasma frequency. 
In these expressions mg and e denote the mass and the charge, respectively. Vap is the alpha-proton drift velocity and 
velocities are expressed in units of the Alfven velocity va = Bq/ iATmmpY/'^. 
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